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Abstract 

The first measurement of diffractive scattering of quasi-real photons with large mo- 
mentum transfer —>■ , where Y is the proton dissociative system, is made using 
the HI detector at HERA. The measurement is performed for initial photon virtualities 
< 0.01 GeV^. Single differential cross sections are measured as a function of W, 
the incident photon-proton centre of mass energy, and t, the square of the four-momentum 
transferred at the proton vertex, in the range 175 < W < 247 GeV and 4 < |t| < 36 GeV^. 
The W dependence is well described by a model based on perturbative QCD using a leading 
logarithmic approximation of the BFKL evolution. The measured |t| dependence is harder 
than that predicted by the model and those observed in exclusive vector meson production. 
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1 Introduction 



The study at the ep collider HERA of exclusive diffractive processes in the presence of a 
hard scale has provided insight into the parton dynamics of the diffractive exchange. The four- 
momentum squared transferred at the proton vertex, t, provides a relevant scale to investigate 
the application of perturbative Quantum Chromodynamics (pQCD) for \t\ ^ ^%^ct> 
this Letter, the first measurement at large t {\t\ > 4 GeV^) of diffractive photon-proton scat- 
tering, 7F, where Y is the proton dissociative system, is presented. The measurement 
is performed at HERA by studying the reaction e^p e^^Y in the photoproduction regime 
with a large rapidity gap between the final state photon and the proton dissociative system Y 
(as illustrated in figure [I^). The centre of mass energy of the system formed by the exchanged 
photon and proton is in the range llh < W < 247 GeV. This process constitutes an extension 
of Deeply Virtual Compton Scattering [2] into the large |t| and small Q"^ regime. 

Diffractive photon scattering can be modelled in the proton rest frame by the fluctuation of 
the incoming photon into a qq pair at a long distance from the proton target. The qq pair is then 
involved in a hard interaction with the proton via the exchange of a colour singlet state. In the 
leading logarithmic approximation (LLA), the colour singlet exchange is modelled by the effec- 
tive exchange of a gluon ladder (figure [IJ5). For sufficiently low values of Bjorken x (i.e. large 
values of W), the BFKL [3] approach is expected to be appropriate for describing the gluon 
ladder. In the LLA BFKL model the gluon ladder couples to a single parton (predominantly a 
gluon) within the proton. The cross section depends therefore linearly on the parton distribu- 
tion in the proton. Due to the quasi-real nature of the incoming photon (Q^ < 0.01 GeV^), the 
transverse momentum of the final state photon, P^, is entirely transferred by the gluon ladder 
to the parton in the proton. The separation in rapidity between the parton scattered by the gluon 
ladder and the final state photon is given by Arj ~ log(s/(Py)^), where s is the invariant mass 
of the system formed by the incoming photon and the struck parton. The proton remnant and 
the struck parton hadronise through fragmentation processes to form the hadronic system Y . 
Assuming parton-hadron duality, hadrons originating from the struck parton correspond to the 
particles with the largest transverse momenta and hence are the closest in rapidity to the scat- 
tered photon. 

The present analysis complements the measurements of exclusive production of p, (f) and 
J /ip mesons at large \t\ [4-7]. The measured W and t dependences of their cross sections were 
found to be in agreement with LLA BFKL based calculations [8-12]. For the process studied 
here, theoretical calculations are simplified by the absence of a vector meson wave function: the 
only non-perturbative part of the calculation is the parton distribution functions of the proton. 
However, the cross section is suppressed relative to that of vector meson production by the 
electromagnetic coupling of the qq pair to the final state photon, making the measurement more 
challenging. 

In the following, measurements of the photon-proton cross sections are presented as a func- 
tion of W and differentially in |t| and are compared to predictions based on LLA BFKL calcu- 
lations [12]. 
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a) 

Figure 1 : a) Schematic illustration of the ep 
process in a LLA BFKL approach. 

2 The HI Detector 




e'yY process, b) Illustration of the 7^*^^ — > 7F 



A detailed description of the HI detector can be found elsewhere [13]. The following briefly 
describes the detector components relevant to this analysis. A right handed coordinate system 
is defined with the origin at the nominal ep interaction vertex, such that the positive z axis 
(forward direction) corresponds to the direction of the outgoing proton beam. The polar angle 
9 and transverse momentum Pt are defined with respect to the z axis. The pseudorapidity is 
defined as 77 = — lntan(^^/2). 

A liquid argon (LAr) calorimeter covers the polar angle range A° < 9 < 154° with full 
azimuthal coverage. The LAr calorimeter consists of both an electromagnetic section and a 
hadronic section. The energy resolution for single particles measured in a test beam isge/E — 
12%/^E/GeV © 1%: for electrons and fT^/E = bQ%/^E/GeV © 2% for hadrons [14]. 
The polar angle region 153° < 9 < 177° is covered by a lead-scintillating fibre calorimeter, 
the SpaCal [15], with both electromagnetic and hadronic sections. The SpaCal has an energy 
resolution for electromagnetic showers oi ge/E = 1%I ^ E jGoN © 1%. 

Charged particles are detected in the polar angle ranges 15° < 9 < 165° by the Central Track 
Detector (CTD) and 5° < ^ < 25° by the Forward Track Detector (FTD). The CTD comprises 
two large cylindrical jet drift chambers, providing precise track measurements in the r — 
plane, supplemented by two 2; -chambers and two multi-wire proportional chambers arranged 
concentrically around the beam-line. The CTD is complemented by a silicon vertex detector 
[16] covering the range 30° < 9 < 150°. The FTD consists of layers of planar and radial drift 
chambers to provide measurement of the 9 and angles of forward tracks, respectively. The 
trackers and the calorimeters are operated within a solenoidal magnetic field of 1.16 T. 

The luminosity is determined from the rate of Bethe-Heitler events, ep — > ep7, measured 
using a Cerenkov crystal calorimeter, the Photon Detector (PD), situated near the HERA beam 
pipes at 2; = —103 m. A second Cerenkov crystal calorimeter, the Electron Tagger (ET), with 
an energy resolution of ge/E = 17%/y/E/GeV © 1%, is located at z = -33 m. The ET 
measures the energy of the positron when scattered through an angle of less than 5 mrad in the 
energy range 8 < E < 20 GeV. The detection of the scattered positron in the ET ensures that 
the exchanged photon is quasi-real with Q"^ < 0.01 GeV^. 
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3 Event Kinematics and Selection 



Following the notation introduced in figure the scattering process, e^p e^'jY, is de- 
scribed by the usual Deep Inelastic Scattering (DIS) kinematic variables: 



p ■ 



where k, p, k' and q are the four-momentum of the incoming positron, the incoming proton, 
the scattered positron and the exchanged photon, respectively. The variable is the virtuality 
of the exchanged photon and y is the inelasticity of the ep interaction, corresponding to the 
relative energy loss of the scattered positron in the proton rest frame. The ep centre of mass 
energy squared is given by s = (k + p^ and the 'jp centre of mass energy squared is W"^ = 
{q + pY — ys- In addition, the longitudinal momentum fraction of the diffractive exchange 
(called the Pomeron in the Regge model [17]) with respect to the proton is defined as: 

q-{p- Py) 

xp = , 

q-p 

where py is the four-momentum of the Y system. The elasticity of the 7p interaction, which can 
be seen as the fractional energy of the exchanged photon transferred to the final state photon, is 
given by 1 — 2/jp, where: 

p-{q- Pj) 

yp = , 

p-q 

Pj being the four-momentum of the final state photon. Finally the square of the four-momentum 
transfer across the diffractive exchange is given by: 

t = {q-p^f = {p-pvf- 



The data used in this analysis were collected with the HI detector during the 1999 — 2000 
running period, when positrons of energy E^, = 27.6 GeV collided with protons of energy Ep = 
920 GeV in the HERA accelerator. The data sample corresponds to an integrated luminosity of 
46.2 pb^^. More details on the present analysis can be found in [18]. 

The data were recorded using a combination of two triggers. Both triggers select an elec- 
tromagnetically interacting particle in the SpaCal which corresponds to the scattered photon 
candidate. One of the triggers requires in addition an energy deposit in the FT, corresponding 
to the scattered positron. The effective trigger efficiency, including time dependent downscale 
factors, amounts to approximately 50%. 

The reconstruction of the kinematic quantities used in the following are approximations 
valid in the limit of small scattering angles of the positron and small transverse momentum 
of the Y system compared to its longitudinal momentum. The quantity y, and hence W, is 
calculated from the scattered positron energy, E^./, measured in the FT using the relation y = 
1 — Ee'/Ee- The relative resolution of W is 4%. To avoid regions of low FT acceptance, the 
energy of the scattered positron is limited to 11 < E^i < 19 GeV, corresponding to 175 < 
W < 247 GeV. In addition, to suppress backgrounds from processes occurring in coincidence 
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with Bethe-Heitler events, it is required that no energy deposits above the noise threshold are 
measured in the PD. 



Photon candidates are selected from energy clusters with small radii detected in the elec- 
tromagnetic section of the SpaCal. If significant energy is measured behind the cluster in the 
hadronic section of the SpaCal, the event is rejected. Events with more than one cluster above 
the noise level in the SpaCal are also rejected. To reduce the background from charged parti- 
cles, the cluster of the photon candidate must have no associated track in the CTD. The photon 
candidates are furthermore required to have an energy Ej > S GeV and a transverse momentum 
P^>2 GeV. 

The hadronic final state Y is reconstructed using a combination of tracking and calorimet- 

ric information. The difference between the total energy E and the longitudinal component of 
the total momentum P^, as calculated from the scattered positron, the final state photon and 
the hadronic system Y, is restricted to 49 < E(E — P^) < 61 GeV. This requirement sup- 
presses non-ep induced background and background due to the overlap during the same bunch 
crossing of a Bethe-Heitler event with a DIS event. For a fully contained ep event the relation 
J:{E- P,) = 2Ee = 55 GeV holds. 

In the case that charged particles from the final state are detected by the tracking detector, 

allowing the primary vertex to be reconstructed, the z coordinate of the vertex has to satisfy 
\z\ < 35 cm. For 25% of the selected events, all charged particles of the final state lie outside 
the detector acceptance and no event vertex is reconstructed. In this case, the time averaged 
vertex position is used for the kinematic reconstruction. 

The kinematic variable |t| is reconstructed as: 

with a relative resolution of 11%. The longitudinal momentum fraction of the diffractive ex- 
change with respect to the proton is reconstructed using the formula: 

where {E + P^)^ is the sum of the energy and longitudinal momentum of the final state photon. 
The event inelasticity of the 7p interaction is reconstructed as: 

- 2{E, - E,,) ' 

where the summation is performed on all detected hadronic final state particles in the event, i.e. 
all measured particles except the scattered positron and the final state photon. This reconstruc- 
tion method has the advantage that the loss of particles along the forward beam pipe only has 
a marginal effect on the reconstruction of yp. Diffractive events are selected by requiring that 
yjp < 0.05. The cut value ensures a large pseudorapidity gap, At^, between the photon and the 
proton dissociative system Y, since ~ e~^'^. 

Finally, to reduce contamination of the signal by non-diffractive background, it is required 
that the difference in pseudorapidity between the photon and the closest final state hadron sat- 
isfies Ar] > 2. This rapidity gap is inferred from the absence of activity in the relevant detector 
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parts, i.e. absence of any track or a cluster of energy deposits above the noise threshold of 400 
MeV in the LAr calorimeter. 

After applying all selection criteria, 240 events remain in the data sample. 



4 Monte Carlo Simulations and Comparison to Data 

Monte Carlo (MC) simulations are used to correct the data for effects of detector acceptance 
and efficiency, to estimate the background and to compare model predictions to the data. All 
generated MC events are passed through the full GEANT [19] based simulation of the HI 
detector and are reconstructed using the same analysis chain as is used for the data. 

The HERWIG 6.4 [20] MC event generator is used to simulate the diffractive high \t \ pho- 
ton scattering using the LLA BFKL model [10-12]. At leading logarithmic accuracy there are 
two independent free parameters in the BFKL calculation: the value of the strong coupling as 
and the scale, c, which defines the leading logarithms in the expansion of the BFKL amplitude, 
ln(iy^/ (c |t| )). In exclusive production of vector mesons, the scale parameter c is related to the 
vector meson mass. In the case of diffractive photon scattering, the unknown scale results in 
the absence of a prediction for the normalisation of the cross section [21]. In the calculations 
considered here [12], the running of as as a function of the scale is ignored. In order to distin- 
guish the as parameter of the LLA BFKL model from the usual strong coupling constant this 
parameter will henceforth be referred to as af^^^. The choice of af^^^ = 0.17 is used for 
the simulation in this Letter. 

In the asymptotic approximation of the calculations [12], the W distribution follows a 
power-law: 

a{W) ~ W^"^", 

where the exponent, also called the Pomeron intercept, is given directly by the choice of af^^'^ 
with ti^o = (Saf^^^/vr) 4 In 2. Using the HERWIG simulation this approximation is found 
to be justified given the current experimental statistical precision. The LLA BFKL model pre- 
dicts an approximate power-law behaviour for the t dependence of the cross section of the form 
d(T/d|t| ~ where n, also predicted by the model, depends only on the parton density func- 
tions (PDFs) of the proton and the value of af^^^. Running coupling effects, not considered 
here, would in principle allow n to depend on t. 

The GRV94 PDFs [22] are used for the proton PDFs in the HERWIG prediction. A com- 
parison with the CTEQ5 [23] and MRST PDFs [24] shows little dependence of the HERWIG 
prediction on the input proton PDFs. 

In order to describe the data, the t dependence of the diffractive photon scattering simulated 
using HERWIG (predicting a value ofn = 3.31 for af^^^ = 0.17) was weighted by a factor 
i.e. the \t\ power is modified from —3.31 to —2.58. This weighted HERWIG prediction 
is used to correct the data for resolution and acceptance effects. 

' The noise level of the LAr calorimeter is about 10 MeV per cell and the average number of cells per cluster 
is typically 60. 
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Possible sources of background are estimated using MC simulations. The background from 
inclusive diffractive photoproduction events {ep eXY, where the two hadronic final states 
are separated by a rapidity gap) is simulated using the PHOJET MC event generator [25]. This 
background contributes when a single electromagnetic particle fakes the photon candidate in the 
SpaCal. It is estimated to amount to 3% of the measured cross section. The background from 
electron pair production (ep ee^e~X) is modelled using the GRAPE event generator [26]. 
This process contributes to the selection if one lepton is detected in the ET, a second lepton fakes 
the photon within the SpaCal and the remaining lepton escapes detection. This background 
contributes 4% of the measured cross section. 

In order to investigate the background from high \t\ diffractive exclusive cu production, 
where the cu decays through the tt+tt^tt^ or 7r°7 channel, a sample was generated using the 
DIFFVM Monte Carlo generator [27]. The contribution from this background process is found 
to be negligible. The background from DIS events, in which the scattered positron fakes the 
photon candidate and an overlapping photoproduction or Bethe-Heitler event gives a positron 
detected in the ET, has been studied and was also found to be negligible. 

In figure[2]the data, corrected for trigger efficiency, are compared to the Monte Carlo simula- 
tions. The background predictions from GRAPE and PHOJET are normalised to the integrated 
luminosity of the data sample. The weighted HERWIG prediction is normalised to the number 
of events obtained after the subtraction of the predicted background from the data. A good 
description of the data by the combined Monte Carlo simulations is observed. 



5 Cross Section Measurement 



The ep e'jY differential cross sections are defined using the formula: 

d'^CTep^e-fY _ N^ata — ^fegr 



where Ndata is the number of observed events corrected for trigger efficiency, N}jgr the expected 
contribution from background events as estimated using the PHOJET and GRAPE Monte 
Carlos simulations, C the integrated luminosity, A the signal acceptance and /SW and At the 
bin widths in W and t, respectively. The acceptance, estimated using the weighted HERWIG 
simulation, is the ratio of the number of events accepted after reconstruction to the number of 
events generated in the defined phase space on hadron level. It accounts for all detector effects, 
including bin-to-bin migrations, as well as geometrical acceptance and detector efficiencies. 
QED radiative correction effects are estimated to be less than 1% [28] and are neglected. 

The 7p — > 7F differential cross section is then extracted from the ep cross section using: 

^f^ = ray)^^ay). (1) 

where the photon flux, r(l^), is integrated over the range Q"^ < 0.01 GeV^ according to the 
modified Weizsacker- Williams approximation [29]. The 7p cross section is obtained by mod- 
elling a^ip-^^Y as a power-law in W, whose parameters are iteratively adjusted to reproduce the 



9 



measured W dependence of the ep cross section. The differential 7p cross section in \t\ is then 
extracted from the ep cross section by correcting for the effect of the photon flux over the visible 
W range (175 < W < 247 GeV). The •yp cross section as a function of W is obtained by first 
integrating equation ([T]) over the \t\ range, and then correcting for the effect of the photon flux 
in each bin in W. More details on the extraction procedure of the 7p cross section can be found 
in [18]. 

The systematic error on the measurement stems from experimental uncertainties and from 
model dependences. They are calculated using the weighted HERWIG simulation of the signal 
process. The sources of systematic error are listed below. For each of them the typical effect on 
the cross section measurement is indicated. 

The experimental systematic errors are: 



• The energy scale uncertainty of ±1% for an electromagnetic cluster measured by the 
SpaCal gives errors in the range of 2% to 4%. 

• The uncertainty of ±2.5 mrad for the measurement of the photon candidate polar angle 
results in an error of up to 3%. 

• The hadronic final state energy scale uncertainty of ±4% leads to an error of less than 
1.5%. 

• The energy scale uncertainty of ±1.5% of the ET produces an error ranging from 1% in 
the highest |t| bin to 10% in the lowest W bin. 

• The uncertainty of ±25% on the noise threshold from the calorimeters gives an error 
varying from 5% at low \t\ to 11% at the highest |t| value. 

• The luminosity is measured with an accuracy of ±1.5% which enters as an overall nor- 
malisation uncertainty. 



The systematic errors due to uncertainty of model parameters are: 



• The uncertainty due to the xp dependence, estimated by weighting the xp distribution 
by the form (1/xjp )^°'^, yields an error varying from 3% in the central \t\ bin up to 9% in 
the lowest W bin. This weight would correspond to changing af^^^ to the values 0.02 
and 0.31. 



• The uncertainty due to the \t\ dependence, estimated by weighting the \t\ distribution by 
the form (l/|t|)='=°-^, leads to an error ranging from 1.5% to 4%. 

• The uncertainty in the modelling of the proton remnant system Y, estimated by weighting 
the My distribution by the form (1/My)^° '^, results in a typical error of 1% to 4%. 

• The uncertainty of 100% assumed on the normalisation of the subtracted inclusive diffrac- 
tive background (from the PHOJET Monte Carlo simulation) leads to an error of approx- 
imately 3% in the highest W bins decreasing to 1% in the highest \t\ bin. 
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• The propagation of the uncertainty on the power-law parameter S in the '-/p cross section 
extraction procedure leads to an error of 4% on d(j^p_^^y/d|t|, independent of and is 
below 1% on a^^^^v- 



The uncertainty on the PHOJET MC normalisation and the model uncertainties from the 
unknown xp , \t\ and My dependences are estimated from data comparisons and defined by 
the range where the model describes the data. Each source of systematic error is varied in the 
weighted HERWIG Monte Carlo within its uncertainty. In each measurement bin, the corre- 
sponding deviation of the normalised cross section from the central value is taken as systematic 
error. The total systematic error is obtained by adding the individual contributions in quadra- 
ture, on a bin-by-bin basis. The largest systematic error on the cross sections comes from the 
uncertainty on the xp and My dependences in the MC simulation. The total systematic error 
on the W dependence of the cross section varies from 10% in the central bins to 17% in the 
lowest bin. The systematic error on the \t\ dependence of the cross section varies from 8% in 
the lowest bin to 14% in the highest bin. An additional global uncertainty of 4% arises from the 
7P cross section extraction procedure. The total systematic errors are comparable to or smaller 
than the statistical errors. 



6 Results 

The cross sections measured for the domain 175 < W < 247 GeV, 4 < |t| < 36 GeV^, 
yiP < 0.05 and Q'^ < 0.01 GeV^ are presented in tabled 

The 7p 7F cross section as a function of W is shown in figure [3l A power-law de- 
pendence of the form a ~ is fitted to the measured cross section. The fit yields 5 = 
2.73 ± 1.02 (stat.)^g 7g (syst.) with x^/n.d.f. = 2.7/2. The contributions from the systematic 
errors are calculated by shifting the data points according to each source of uncertainty, taking 
correlations into account, and repeating the fit. The errors are then added in quadrature to obtain 
the total systematic error. 

The steep rise of the cross section with W is usually interpreted as an indication of the pres- 
ence of a hard sub-process in the diffractive interaction and of the applicability of perturbative 
QCD. The present 6 value, measured at an average |t | value of 6. 1 GeV^, is compatible with that 
measured by HI in diffractive J/ip photoproduction of 5 = 1.29 ± 0.23(stat.) ± 0.16(syst.) [6] 
at an average |t| of 6.9 GeV^. The Pomeron intercepts associated with these 5 values correspond 
to the strongest energy dependences measured in diffractive processes. 

The 7p cross section differential in \t\, at W = 219 GeV, is shown in figure HI Figure |4] 
also shows the comparison of the cross section to a fit of the form da/dt ~ |t| The fit result 
is n = 2.60 ± 0.19 (stat.)^QQg (syst.) with x^/n.d.f. = 1.6/1. The |t| dependence is harder 
than that measured by HI in the diffractive photoproduction of J/^ mesons at large |t| [6] 
corresponding to n = 3.78 ± 0.17(stat.) ± 0.06(syst.). 

In figures[3]and|4]the measured cross sections are compared to predictions of the LLA BFKL 
model, using the HERWIG Monte Carlo, as described in sectionlH with no \t\ weighting. The 
predictions are normalised to the integrated measured cross section, as the normalisation is 
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HI measured ep e'jY cross sections 


W 
[GeV] 


W range 
[GeV] 


dcTep^e-yY / dW 

[pb/GeV] 


T{W) 
[GeV-^] 


(^■yp^'yY 

[nb] 


185 


175 - 193 


0.414 ±0.069 ±0.072 


0.0565 


2.02 ±0.34 ±0.35 


202 


193 - 211 


0.318 ±0.046 ±0.033 


0.0431 


1.86 ±0.27 ±0.19 


220 


211 - 229 


0.434 ±0.062 ±0.051 


0.0328 


3.06 ±0.44 ±0.36 


240 


229 - 247 


0.404 ± 0.080 ± 0.044 


0.0246 


3.48 ± 0.69 ± 0.38 


\t\ 
[GeV2] 


\t \ range 
[GeV2] 


dcTep^e-yv/dltl 

[pb/GeV2] 


T{W = 219 GeV) 
[GeV-i] 


dcr^p^^Y /d\t\ 
[pb/GeV2] 


6 


4.0-8.3 


4.04 ±0.42 ±0.36 


0.0333 


401 ± 41 ± 36 


12 


8.3-17.3 


0.58 ±0.08 ±0.06 


0.0333 


57.8 ±8.1 ±6.2 


25 


17.3-36.0 


0.13 ±0.03 ±0.02 


0.0333 


12.5 ±3.1 ±1.8 



Table 1: The cross sections for the processes ep — > e^yY and the 7p 7F, measured in the 
range yp < 0.05 and Q"^ < 0.01 GeV^. The upper part of the table presents the measured cross 
sections for different values of W at an average {\t\) = 6.1 GeV^. The lower table presents 
the measured cross sections differential in \t\ atW = 219 GeV The first errors are statistical, 
the second systematic. The photon flux T and corresponding ranges in W and \t\ used for the 
measurements are also given. 



not predicted by the LLA BFKL calculation [21]. The W dependence of the cross section 
is well described by the LLA BFKL prediction, as shown in figure [3l The sensitivity of the 
measurement to the free parameter af^^^ in the theoretical prediction is illustrated in figures |3] 
and m Using 5 = Auq = A (3af ' */ vr) 4 In 2, the measured W dependence leads to af '* = 
0.26 ± 0.10 (stat.)+J;j7 (syst.). Predictions are shown in figureHfor the values af = 0.14 
and 0.37. The LLA BFKL curve corresponding to af^^^ = 0.26 coincides with the solid fit 
line. 

Previous measurements of diffractive scattering at HERA are well described by BFKL 
predictions with af^^^ values similar to the value measured in this Letter. ZEUS measure- 
ments of exclusive p, and J/'ip production at large |t| are best described with the value of 
af^^^ = 0.20 [4]. The HI measurement of high \t\ p production [5] is compatible with 
^BFKL _ g 2Q ^nd the HI measurement of high \t\ J/ if) production [6] is described using 
Q,BFKL _ Q ig_ j]^e LLA BFKL prediction with af^^^ = 0.17 gives a good description of the 
double dissociation process with a rapidity gap between jets measured by HI [30]. Events with 
rapidity gaps between jets were also measured by ZEUS [31] and are found to be compatible 
with a model which uses af^^^ = 0.11. Note, however, that for these measurements the hard 
scale, corresponding to the jet transverse momentum squared, is of the order of or larger than 
20 GeV^. 

As shown in figure H the LLA BFKL calculation for af^^^ = 0.14, 0.26 and 0.37, aU of 
which give a reasonable description of the W dependence, predict steeper |t| distributions than 
is measured in the data. The same effect cannot be established for the exclusive p measurement 
[5], where the measured t range is limited to |t| < 8 GeV^, although an underestimate of the 
cross section was observed at the largest values of The present situation is in contrast with 
the analysis of J/ip production [4,6], where the \t\ dependence was found to be well described 
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by the LLA BFKL prediction over a similar range in t. 



7 Conclusions 

Using the HI detector, diffractive photon scattering, 7p — > 7F, where the final state photon 
carries a large transverse momentum and is well separated from the proton dissociative sys- 
tem Y, is measured for the first time at HERA. The measurement provides a unique test of the 
underlying QCD dynamics of the diffractive exchange. 

Cross sections are presented as a function of W and differentially in |t | . A fit of the form 
performed on the cross section gives 6 = 2.73 ± 1.02 (stat.)^Q7g (syst.). This strong energy 
dependence is compatible with that measured for the exclusive diffractive J/^ production at 
high \t\. A fit of the form yields n = 2.60 ± 0.19 (stat.)l'iQ'Qg (syst.), corresponding to a 
harder |t| dependence of the cross section than measured for high |t| diffractive J/^ production. 

The measured cross sections are compared to the predictions of an LLA BFKL model. A 
good description of the W dependence of the cross section is found, whereas the LLA BFKL 
model predicts a\t\ dependence that is too soft and hence unable to describe the data. 
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Figure 2: Distributions of the selected events as a function of a) the final state photon energy, 
b) the final state photon polar angle, c) the scattered positron energy, d) e) xjp , and f) yjp. 
The data corrected for trigger efficiency (black points) are compared with the simulation of 
diffractive high \t\ photons from the weighted HERWIG based on the LLA BFKL calculation 
(open histogram), and the background from inclusive diffractive photoproduction simulated 
with PHOJET (full histogram) and dilepton production simulated with GRAPE (hatched his- 
togram). The HERWIG prediction is normalised to the number of data events after background 
subtraction. The total systematic uncertainty of the simulation is shown by the dark grey shaded 
band. 
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Figure 3: The 7p cross section of diffractive scattering of photons as a function of W in the 
phase space defined by 4 < |t| < 36 GcV^ , yjp < 0.05 and < 0.01 GcV^. The average 
value is {\t\) = 6.1 GeV^. The inner error bars show the statistical errors and the outer error 
bars show the statistical and systematic errors added in quadrature. The solid line shows the 
result of a fit to the cross section of the form with d = 2.73 ± 1.02 (stat.)^[J'^g (syst.). This 
line also corresponds to the LLA BFKL model prediction from the HERWIG event generator 
with af^^^ = 0.26. Two additional curves show the LLA BFKL predictions for the additional 
choices of af^^^ — 0.14 and 0.37 corresponding to one standard deviation from the fit value. 
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Figure 4: The 7p cross section of diffractive scattering of photons differential in \t\ for W = 
219 GeV, yjp < 0.05 and < 0.01 GeV^. The inner error bars show the statistical errors 
and the outer error bars show the statistical and systematic errors added in quadrature. The 
solid line shows the result of a fit to the cross sections of the form with n = 2.60 ± 
0.19 (stat.)^g Qg (syst.). Three additional curves show the LLA BFKL model predictions from 
the HERWIG event generator for the values af^^^ = 0.14, 0.26 and 0.37. 
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